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We conjecture that the confinement- deconfinement phase transition in QCD at large number 
of colors TV and Nf <C N at T 7^ and /i 7^ is triggered by the drastic change in 9 behavior. 
The conjecture is motivated by the holographic model of QCD where confinement -deconfinement 
phase transition indeed happens precisely at the value of temperature T — T c where 9 dependence 
experiences a sudden change in behaviorJJ]. The conjecture is also supported by quantum field theory 
arguments when the instanton calculations (which trigger the 9 dependence) are under complete 
theoretical control for T > T c , suddenly break down immediately below T < T c with sharp changes 
in the 9 dependence. Finally, the conjecture is supported by a number of numerical lattice results. 

qq ■ We employ this conjecture to study confinement -deconfinement phase transition of dense QCD at 

large fi in large N limit by analyzing the 9 dependence. We find that the confinement- deconfinement 
phase transition at Nf <C TV happens at very large quark chemical potential fx c ~ VNAqcd- This 

£S| ■ result agrees with recent findings by McLerran and Pisarski[2(|. We also speculate on case when 

' Nt ~ N. 
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. I. INTRODUCTION 

i-C ' Understanding the phase diagram at nonzero external parameters T, fi is one of the most difficult problem in QCD. 

Obviously, this area is a prerogative of numerical lattice computations. However, some insights about the basic 
O ■ features of the phase diagram may be inferred by using some analytical approaches. In particular, some qualitative 
questions can be formulated and answered by considering a theory with large number of colors N or/and flavors Nf, 
see recent papers on then subject [J 0, 9 and references on previous works therein. Generically, to study a phase 
transition one should find an appropriate order parameter. It is easy to find an order parameter for gluodynamics 
£SJ . when light quarks are not present in the system. If massless quarks are introduced into the system, one can study a 
chiral phase transition and use the chiral condensate as an order parameter. For massive, but light quarks this is not 
an option. However, in the limit of very large TV one can consider the free energy as an order parameter. In confined 
phase it is order of one, while in deconfined phase it is order of ~ TV 2 . Small number of flavors Nf -C TV (massless or 
massive quarks) does not change the basic picture. 

We formulate a different criteria for confinement -deconfinement phase transition, and therefore we use a different 
order parameter to analyze the phase transition. The new criteria is based on observation that the deconfined phase 
transition is always accompanied by very sharp changes in 9 behavior which represents our basic conjecture. Therefore, 
in principle, if our conjecture is correct, one can use any order parameter which nontrivially depends on 9 and study 
this dependence on two sides of the phase transition line. Very natural question immediately comes into mind: 
why and how these two different things (phase transition vs sharp 9 changes ) could be linked? What is the basic 
motivation for this proposal? First of all, this criteria is motivated by the observation that in holographic model of 
QCD the confinement -deconfinement phase transition happens precisely at the value of temperature T = T c where 9 
dependence experiences a sudden change in behavior Secondly, the proposal is supported by the numerical lattice 
results [3] -0, see also a review article [§], which unambiguously suggest that the topological fluctuations are strongly 
suppressed in deconfined phase, and this suppression becomes more severe with increasing TV. These general features 
observed in the lattice simulations have very simple explanation within our proposal on the origin of the confinement 
-deconfinement phase transition, see next section for details. Finally, our new criteria is based on a physical picture 
which can be shortly summarized as follows. 

For sufficiently high temperatures T > T c the instanton gas is dilute with density ~ e -i( T ) N which implies a strong 
suppression 1 of the topological fluctuations at large TV where 7(T) > 0, see below for details on structure of j(T)— 
function. The calculations in this region are under complete field theoretical control and the vacuum energy has a nice 
analytic behavior ~ cos9e~~ ,l - T ' ,N as function of 9. At the critical value of temperature, T = T C where j(T) changes 
the sign, the instanton expansion breaks down and one should naturally expect that at T = T c there should be a sharp 
transition in 9 behavior as simple formula ~ cos 9 can only be valid when the instanton gas is dilute and semiclassical 



See [loll and references therein for earlier discussions on the subject. 
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calculations are justified which is obviously not the case for T < T c . Therefore, it is naturally to associate sharp 
changes in 9 behavior with confincment-deconfincmcnt transition, just as in the holographic model[l]. There is a very 
narrow window of temperatures in deconfmed phase, < (T — T c )/T c < l/N when the instanton expansion is not 
valid. This vicinity of T c is extremely interesting, see our comments about physics in this region in conclusion. This 
region shrinks to a point at N = oo. 

The main goal of this paper is to apply this criteria to the region with large chemical potential at large N and 
Nf <C N and make a specific prediction on magnitude fJ> c (T) for confinement-deconfinement transition line at large 
\x and sufficiently small T <C /x. The corresponding estimation of /x c (T) is based on well-developed instanton calculus 
in deconfmed phase where dilute gas approximation is justified. 

The plan of the paper is as follows. We start in Section II by reviewing recent work [l[ on estimation T c using 
instanton calculus. We also present a picture explaining how and why two apparently different phenomena (sharp 
changes in 9 and confincment-deconfincmcnt transition) may in fact be tightly linked. In section III we apply the 
same technique to argue that the confinement- deconfinement phase transition happens at very large quark chemical 
potential [i c ~ ^/NAqcd, where fi — [ib/N is already properly scaled quark chemical potential. This result agrees 
with recent analysis by McLerran and Pisarski Q which was based on fundamentally different starting point. Finally, 
in section IV we make few comments for the case when number of flavors Nf ~ N. 



II. CONFINEMENT- DECONFINEMENT PHASE TRANSITION IN HOT QCD AT LARGE N. 

We start with a short review of ref.[l[ where the conjecture (that the confinement-deconfinement phase transition 
happens precisely where 9 behavior sharply changes) was implemented for large N QCD at T ^ 0. Such a sharp 
transition is indeed observed in the holographic model of QCD. From quantum field theory viewpoint such a transition 
can be understood as follows. Instanton calculations are under complete theoretical control in the region T > T c as 
the instanton density is parametrically suppressed at large N in deconfmed region[l|, 

W0)~e-^cos0, 7 = \±±l n (-^-) -1.861. (1) 

L -J \ A QC-D/ J 

It is assumed that a higher order corrections may change the numerical coefficients in j(T), but they do not change the 
structure of eq. ([1]) . The critical temperature is determined by condition 7 = where exponentially small expansion 
parameter e~~ lN suddenly blows up and becomes exponentially large. Numerically, it happens at 
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= => T C (N = oo) ~ 0.53Aqcd, (2) 



where Aqcd is defined in the Pauli -Villars scheme. Our computations are carried out in the regime where the 
instanton density ~ exp(— 7 AT) is parametrically suppressed at any small but finite 7(T) = e > when N = 00. From 
eq. {]]) one can obtain the following expression for instanton density in vicinity of T > T c , 

V inst (9) ~ cosfl ■ e- aN (^), 1 » (j^r^j » l/N. (3) 

where a = and T C (N = 00) ~ 0.53Aqcd are estimated at one loop level. Such a behavior does imply that the 
dilute gas approximation is justified even in close vicinity of T c as long as T ~ Tc ^> Therefore, the 9 dependence, 
which is sensitive to the topological fluctuations is determined by ((3|) all the way down to the temperatures very close 
to the phase transition point from above, T = T c + 0(1/N). The topological susceptibility is order of one for T < T c 
in confined phase while it vanishes ~ e~ 7jv — > for T > T c in deconfmed phase. Non topological quantum fluctuations 
on the other hand could be quite large in this region, but they do not effect the structure of eq. ([3]). We do anticipate, 
of course, that the perturbative corrections in the instanton background may change our numerical estimate for T c 
and a. However, we do not expect that a qualitative picture of the phase transition may be affected as a result of 
these corrections. We note that the lattice numerical computations [|[ -0] do suggest that the topological fluctuations 
are strongly suppressed in deconfmed phase immediately above T c , and this suppression becomes more severe with 
increasing N starting from physically relevant case N = 3. Holographic QCD also supports this picture[l|. We do 
not expect any changes in the picture when small number of flavors Nf -C N are introduced into the system 2 . 



2 We have to make the following remark here in order to avoid any confusions later in the text. In the presence of the massless chiral 
fermions the 8 dependence goes away in QCD in both phases: confined as well as deconfined. It is a simple reflection of the fact that 
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There are three basic reasons for a generic structure (11I2I3P to emerge: 

1. The presence of the exponentially large "T— independent" contribution ( e.g. e + 1 -^6N j n (JTJ)) . This term 
basically describes the entropy of the configuration. It is due to a number of contributions such as a number of 
embedding SU(2) into SU(N) etc; 

2. The presence of the "T— dependent" contribution to Vi ns t(0) which comes from J n(p)dp integration, see below 
([8|). It is proportional to 



A 



QCD 



-N 



= cxp 



11„ , / nT 
N ■ In 
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(4) 



3. The fcrmion related contributions such as a chiral condensate, diquark condensate or non-vanishing mass term 
enter the instanton density as follows ~ ({pip) N t ~ e JV '( Kln IW'V>}|) i p or K = Ej_ _> o this term obviously leads to a sub 
leading effects 1/N in comparison with two main terms in the exponent ([1]). Therefore, such terms can be neglected 
as they do not change any estimates at N — oo. It is in accordance with the general arguments suggesting that the 
fundamental fermions can not change the dynamics of the relevant gluon configurations as long as Nf -C N. 

The crucial element in this analysis is that both leading contributions (items 1 and 2 above) have exponential e N 
dependence, and therefore at N — > oo for T > T c the instanton gas is dilute with density e~ lN , 7 > which ensures 
a nice cos 9 dependence ([3]) , while for T < T c the expansion breaks down, and 9 dependence must sharply change at 
T <T C . We have identified such sharp changes with first order phase transition. 

Once T c is fixed one can compute the entire line of the phase transition T c (/i) for relatively small /1 <C T c for large 
but finite N ^> Nf. The result in the leading loop order can be presented as follows 



T c (/x) = T c (/z = 0) 



4iV7r 2 T c 2 (/j = 0) 



H < ttT c , Nf < N. 



(5) 



As expected, fi dependence goes away in large N limit in agreement with general large N arguments 11| . This 
formula is in excellent agreement with numerical computations (T2I Il3l . [l4| which show very little changes of the 
critical temperature T c with \i for sufficiently small chemical potential. In particular, even for the case Nf = 2, N = 3 
where the expression ([5J is not expected to give a good numerical estimate, it still works amazingly well even for 
N = 3. Indeed, the result quoted in [12j | can be written as 



T c (nY at = Tjji = 0)' Qf 1 - 0.500(67) 



TT 2 T?(p = 0) lat 

It should be compared with our theoretical prediction (J5]) for this case 



N f = 2, N = 3. 



T c {n) th = T c (p = 0) 
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The eq.(JS|) suggests very slow change of T c with /i at large N. Such slow variation implies that a sufficiently large 
changes of order one AT C ~ ^qcd may occur only when chemical potential changes are very large, A/i ~ s/N Aqcd ■ In 
next section we confirm this expectation by a direct computations of /x c (T = 0) where we predict that the confinement 
-deconfinement phase transition happens at very large /z c (T = 0) ~ \/NKqcd if Nf <C N. 

One more comment on this proposal. Our conjecture (that the confinement- deconfinement phase transition in QCD 
is triggered by the drastic change in 9 at the same point T = T c ) implicitly implies that the configurations which are 
responsible for sharp 9 changes must also play a significant role in confined phase at T < T c . On the other hand, at 
T > T c the dilute instantons completely determine the 9 dependence ([3|) while at T < T c the small size instantons 
obviously can not provide confinement [15j. How can this be consistent with our conjecture that these two things must 



one can redefine the fermi fields in the chiral limit such that 8 parameter completely disappears from the partition function. To avoid 
the identical vanishing of Vi ns t{9) one can introduce a non-zero quark mass m q ^ 0. It does not effect any of our estimates as long as 
Nf -C N as all such changes lead to a sub leading 1/N corrections, see item 3 below. Our goal here is to study the coefficient in front 
of cos 8 in deconfined phase. By such an analysis we trigger the point when this coefficient is suddenly blows up, and the 8 dependence 
must drastically change. The sharp changes of this coefficient ~ Vi ns t (8) we identify with complete reconstruction of the ground state, 
drastic changes of the relevant gluon configurations, and finally, with confinement- deconfinement phase transition. One should also 
remark here that the assumption made in [l| on non-vanishing chiral condensate in vicinity T > T c as a holographic model of QCD 
suggests, is not crucial for our arguments to hold as it leads to a sub leading 1/7V correction, see item 3 below. 
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be linked? We note that quark confinement can not be described in the dilute gas approximation, when the instantons 
and anti-instantons are well separated and maintain their individual properties (sizes, positions, orientations), as it 
happens at large T > T c . However, in strongly coupled theories the instantons and anti-instantons lose their individual 
properties (instantons will "dissociate ") their sizes become very large and they overlap. The relevant description is 
that of instanton-quarks 3 , the quantum objects with fractional topological charges ztl/N which become the dominant 
quasi-particles. The instanton quarks carry, along with fractional topological charges, the fractional 1/N magnetic 
charges which are capable to propagate far away from instantons- parents being strongly correlated with each other. 
For such configurations the confinement is a possible outcome of the dynamics. It makes the instanton quarks to 
become the perfect candidates to serve as the dynamical magnetic monopoles, the crucial element of the standard 
't Hooft and Mandelstam picture for the confinement [HI, . This basically represents our proposal for the answer 
on the question formulated above. One should emphasize that our arguments (that the instantons dissociate into 
instanton quarks in confined phase) are not based on any semiclassical analysis performed in strongly coupled regime. 
Rather, they are based on the following observation. The dimensionality of the moduli space of the relevant statistical 
ensemble precisely coincides with the corresponding k instanton measure 4Nk with k being integer. This "coincidence" 
holds for any gauge group G, not limited to SU(N) case (2£j. Similar picture on dissociation of the instantons into 
instanton quarks in confined phase has been also recently advocated in [2lj . [22| . see our comments on these papers 
inQ. 

The main lesson of this section can be formulated as follows: we presented a number of arguments suggesting 
that our proposal (which relates two naively unrelated phenomena: confinement- deconfinement phase transition and 
drastic changes in 9 behavior) is consistent with all previously known studies. In particular, it includes: 

• lattice computations of topological susceptibility in the vicinity of the phase transiti on^ , | a , Jfj. El § ; 

• lattice computations of the critical temperature T c (/j,) as a function of ^ at small /i [12. Il3l. Il4j|; 

• analysis of the holographic models of QCD in vicinity of the phase transition at nonzero temperature [l], H^, [24[ ■ 
In such circumstances when the outcome which follows from the basic conjectured principle agrees with all known 

results, one should naturally try to extend the corresponding analysis to the regions in the parametrical space which 
are presently not accessible for study by other means. To be more specific, we want to analyze the confincmcnt- 
deconfinement phase transition at very large chemical potential fj, > Aqcd when available technique does not allow 
to perform the lattice computations. One should also note that presently available holographic models of QCD also 
can not address this question. With this motivation in mind we want to analyze confinement- deconfinement phase 
transition at large chemical potential and compare the obtained results with corresponding analysis [2|, |3J] which is 
based on fundamentally different starting point. 



III. CONFINEMENT- DECONFINEMENT PHASE TRANSITION IN DENSE QCD AT LARGE N. 

In this section we estimate the value of [i c where the instanton expansion breaks down and therefore, the 9 depen- 
dence should experiences a sharp change. According to our conjecture we should identify this place with the phase 
transition point. Similar arguments have been put forward previously [25] for numerical estimation of /i c for small 
N,Nf = 2,3, see also review talk on this subject 26]. Our goal here is quite different: we want to understand an 
analytical dependence of /j, c (N,Nf) as a function of N,Nf at very large N and finite Nf <C N in order to compare 
with results of refs. @, Q where the authors presented a very strong argument suggesting a very large /Lt c ~ \/N where 
the phase transition could happen. 

We follow the same logic as inQ, and study the 9 dependence in order to make a prediction about the phase transition 
point /i c . In the regime ji > /i c the 9 dependence is determined by the dilute instanton gas approximation. We expect 
that the expansion breaks down only in close vicinity of [l c at large N as it happens in our previous analysis with 
phase transition at T = T c . According to the conjecture this point will be identified with confinement- deconfinement 
phase transition point fj, c . In the present case of analyzing fj, c rather than T c discussed previously [l(, we do not 
have any support from the lattice computations, nor from holographic models. Still, the basic governing principle 
remains the same. Therefore we identify the point where instanton expansion breaks down (and correspondingly a 
point where a simple cos 9 sharply changes to something else) with the point /i c where the phase transition happens. 



Instanton quarks originally appeared in 2d models. Namely, using an exact accounting and resummation of the n-instanton solutions 
in 2d C'P 1 ^^ 1 models, the original problem of a statistical instanton ensemble was mapped unto a 2d-Coulomb Gas (CG) system of 
pseudo-particles with fractional topological charges ~ 1/jV |16H . This picture leads to the elegant explanation of the confinement phase 
and other important properties of the 2d CP N—1 models Il6ll . Unfortunately, similar calculations in id gauge theories is proven to be 
much more difficult to carry out [l7| . 
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In our estimates below we assume that the color superconducting phase is realized in deconfined phase for all N, see 
e.g. recent review |27j . It is known though that for extremely large N — oo one could expect that another pha se is 
more energetically favorable 28]. Still, for all reasonably large N the color superconducting phase prevails |29l|. In 
any case, the difference between the two options would lead to a sub- leading 1/N corrections as explained in item 3 
above. 

As we shall see below, the instanton density in deconfined phase has the following generic behavior, ~ 
cos 8 exp [— Nj ([i)], where 7(/i) ~ const. + 0(1/ N) in large N limit. Such a behavior implies that for any small 
(but finite) positive 7 > the instanton density is exponentially suppressed and our calculations are under complete 
theoretical control. In contrast: at arbitrary small and negative 7 < the instanton expansion obviously breaks down, 
theoretical control is lost as an exponential growth ~ exp(|7|iV) for the instanton density makes no sense. The 8 
behavior must drastically change at this point. Therefore, the value of p c is determined by the following condition, 

7(M = Mc)=0 => p c = cA QC D- (6) 

Our goal is to compute the coefficient c by approaching the critical point p c from deconfined side of phase boundary. 
Therefore, we will be interested in the instanton density in the dilute gas regime at p > p c where analytical instanton 
calculations are under control. 

As we already mentioned in footnote 2, pg.2 the 8 dependence goes away in full QCD in both phases: confined 
as well as deconfined in the presence of the massless chiral fermions. However we are interested in the magnitude 
of the instanton contribution ~ Vmst(Q) in deconfined phase rather than in 8 dependence of full QCD. Precisely this 
coefficient triggers the point where the instanton expansion suddenly blows up. The sharp changes in Vi ns t(8) we 
identify with complete reconstruction of the ground state, drastic changes of the relevant gluon configurations, and 
finally, with confinement- deconfinement phase transition. To avoid identical vanishing of Vinst(^) m the presence of 
massless fermions one can assume a non zero chiral condensate in deconfined phase, as it has been done for hot matter 
in[l[ with motivation from holographic model of QCD. Also, one can assume a non-vanishing masses m q 7^ for the 
fermions, or non- vanishing diquark condensate (tpt/>) ^ to avoid identical vanishing of Mnst(^) for dense matter at 
large p. None of these assumptions effects any numerical estimates given below in the limit N — > 00, Nf <C N, as all 
these assumptions lead to a sub-leading effects ~ 1/N which will be ignored in what follows. We shall see in section 
IV that such kind of assumptions indeed play a crucial role but only when Nf ~ N. 

To be definite, we assume that the non-vanishing diquark condensate (iptf?) ^ develops for a > p c . A precise 
magnitude of the diquark condensate is not essential for our calculations as it effects only sub- leading terms ~ 1 jN 
which will be consistently ignored in what follows. The instanton-induced effective action for Nf massless fermions 
can be easily constructed. In particular, for Nf = 2 flavors, u, d the corresponding expression takes the following 



form, [3j,l3J,l3J,l3J,l3J,l35|, 



L lu>l =<-' ! 'j<i,>n (/,)[ ^ttV) {U,,«>l)UIi,«Il)+ !' ! 



_3_ 

32 



3 

(u R X a u L )(dRX a d L ) - ^{uR<J fl v\ a u L )(d R <j t2l/ \ a d L ) 



H.c. 



We wish to study this problem at nonzero chemical potential u and nonzero small temperature T <C p, (to be discussed 
later in the text). We use the standard formula for the instanton density at two- loop order (30l. l3ll. [32l. l33l. l3~il |3"H 

n(p) = C N (f3 I (p)) 2N p- 5 exp[-f3 II (p)} x exp[-(N f p 2 + ±(2N + N f )n 2 T 2 )p 2 }, (8) 



where 



0.466e- 1679iV 1.34 JV / a . . , , . . . Q . . „ . . b> , (2^(p) 

Cn= r N _ 1 yj N _ 2 )\ ' ^ = -^MMqcd), Mp) = frO>) + ^ log 



3 3 ; 3 3 N 



This formula contains, of course, the standard instanton classical action exp(— 8ir 2 /g 2 (p)) ~ exp[— Pi(p)] which 
however is hidden as it is expressed in terms of Aqcd rather than in terms of coupling constant g 2 (p) . The chemical 
potential p = pb/N in this expression is already properly normalized quark chemical potential (rather than baryon 
chemical potential). By taking the average of e q.([71) over the state with nonzero vacuum expectation value for the 
diquark condensate (ipip) ^ as described in [25l |36|. integrating over p, and taking large N limit using the standard 
Stirling formula 

1 1 



T(N + 1) = V2^NN N e- N ( 1 + — + 0(— ) ) (9) 



6 



one finds the following expression for the instanton induced potential 4 , 




V inst (9) ~e-~<"cos9, 7= [— In I I -l.lj, ^ = Nfi z , (10) 

where we introduced reduced chemical potential p, = fi/yN and neglected all powers N p in front of e~~' yN . The crucial 
difference in comparison with similar computation at nonzero temperature {T]) is emerging of parameter fl instead 
of the original quark chemical potential fi = \f~N \x. It implies that the critical chemical potential where 7 changes 
the sign (and therefore where the phase transition is expected) is parametrically large fi c ~ y/N because Ji c ~ 1, 
see below for numerical estimates. The origin for this phenomenon can be traced from eq. ((8]) where temperature 
dependent factor in the instanton density is proportional to ~ N while chemical potential enters this expression with 
factor ~ Nf <C N. Therefore, a very large chemical potential /1 ~ \f~N Aqcd is required in order to achieve the same 
effect as temperature T ~ Aqcd- The physics of this phenomenon can be explained as follows: at T ~ 1 a large 
number of gluons ~ A^ 2 can get excited while at [X ~ 1 only a relatively small number of quarks in fundamental 
representation ~ N can get excited. Therefore, it requires a very large chemical potential fi 2 ~ TV in order for 
fundamental quarks play the same role as gluons do at T ~ 1. As explained above, the critical chemical potential is 
determined by condition 7 = where exponentially small expansion parameter e~' yN at /i > fj, c suddenly blows up at 
fi < /i c . Numerically, it happens at 

- 1.1 





fi c (N = 00) ~lA-A Q CD X Tr , N f <N, (11) 



where Aqcd is defined in the Pauli -Villars scheme. The topological susceptibility vanishes ~ e~ lN — > for fj, > fi c 
while it must be drastically different for /i < fj, c as 9 dependence must experience some drastic changes in this region 
as the instanton expansion breaks down, and therefore simple cos 9 dependence must be replaced by something else. 
It is very likely that the standard Witten's arguments (valid for the confined phase) still hold in this region /1 < fi c in 
which case the topological susceptibility is order of one. 

The Aqcd in the Pauli -Villars scheme which enters our formula (fTTj) is not well-known numerically. Therefore, for 
numerical estimates one can trade Aqcd m favor of T C (N = 00) at \i = estimated in[l|, see eq. @. Therefore, our 
final numerical estimate for fj, c (N — 00) can be presented as follows, 




fi c (N = 00) ~ 2.6 • i / —— ■ T C (N = 00,11 = 0), Nf N. (12) 



If one uses the numerical value for T C (N = 3) ~ 260 MeV [5, 8], one arrives to ^ C (N = 00) ~ 690y/N/N f MeV which 
is our final numerical estimate for the critical chemical potential where deconfined phase transition is predicted for 
very large N . Few remarks are in order: 

a. The most important result of the present studies is the observation that the confinement- deconfinemcnt phase 
transition according to (fTTj) happens at very large [i c ~ \[N if Nf <C N . This is consistent with the results of Q 
where parametrically large scale for /i c ~ \[N had been predicted. However, the technique of ref. does not allow to 
answer the question whether the transition would be the first order or it would be a crossover. Within our framework 
at ^> 1 and Nf <C A^ the entire phase transition line (which starts at T = T c ~ Aqcd at fi = and ends at 
(j, = fi c ~ ^/NAqcd at T = 0) is predicted to be the first order phase transition at large N and Nf -C N. This 
is because the nature for the phase transition along the entire line is one and the same: it is drastic changes of 9 
dependence when the phase transition line is crossed. 

b. Our computations are carried out in the regime where the instanton density ~ exp(— 7 AT) is parametrically 
suppressed at N = 00. From eq. (|10[) one can obtain the following expression for instanton density in vicinity of 

M > Mo 

W0)~cos0.e-^(^), i« «l, (13) 



4 The diquark condensate in large N limit has behavior (tp^) ~ ex P( — ~) ~ cxp(— V^V), see e.g. review [27|1 . It is still a sub-leading 

1/y/N effect in comparison with the main terms II10L Author thanks an anonymous referee for pointing out on this, potentially large, 
correction. 
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where a is 11/3 at one loop level, but the perturbative corrections could be large and they may considerably change 
this numerical coefficient. Such a behavior (|13[) does imply that the dilute gas approximation is justified even in 
close vicinity of p c as long as M ~ Mc 3> jf. In this case the diluteness parameter remains small. We can not rule out, 
of course, the possibility that the perturbative corrections may change our numerical estimate for p c . However, we 
expect that a qualitative picture of the phase transition advocated in this paper remains unaffected as a result of 
these perturbative corrections in dilute gas regime. 

c. In our estimate for p c we neglected (log pAqcd) 1 * in evaluating of the j dp integral. The corresponding correction 
changes our estimate (|lip very slightly, and it will be ignored in what follows. Numerical smallness of correction is 
due to the strong cancellation between the second loop contribution in the exponent (term proportional to V /b) and 
the first loop contribution in the pre-exponent in eq. ([8]). 

d. Once p c is fixed one can compute the entire segment of the phase transition line p c {T) for relatively small T. 
Indeed, in the dilute gas regime at p > p c the T dependence of the instanton density is determined by a simple 
insertion ~ exp[-2/3Nn 2 T 2 p 2 ] in the expression for the density |8]), In the leading loop order p c (T) varies as follows, 

r Nn 2 T 2 i 

Mc (T)= Mc (T = 0)[l- 3JV ^ (T = 0) ], VNT«, C . (14) 

One should remark that a variation of the critical chemical potential Ap c (T) is very large ~ y/N when the temperature 
variation AT ~ 1 is order of one in units of Aqcd- This is in huge contrast with a similar expression |5]| which shows 
very little change ~ 1/N of the critical temperature AT c (/i) ~ X/N with variation of chemical potential of order one, 
A/i ~ 1. The nature of this difference between p c and T c was already mentioned before and can be explained by the 
fact that at T ~ 1 a large number of gluons ~ TV 2 can get excited while at p, ~ 1 only a relatively small number 
of quarks in fundamental representation ~ N can get excited. Therefore, it requires a very large chemical potential 
p? ~ N when quarks can play the same role as gluons do at T ~ 1 as long as Nf <C N. 



IV. DECONFINEMENT TRANSITION IN HOT AND DENSE QCD AT N f ~ N. SPECULATIONS. 

Our estimations (|11|12I14|) have been derived under assumption that N — > oo while Nf is kept fixed such that 
k = Nf/N — > 0. In particular, for the case of hot matter with T ^ studied in [l[ the fermi fields and the chiral 
condensate were introduced exclusively with a single purpose to elucidate the physical interpretation of the phase 
transition for pure gluodynamics rather than for full QCD. The physical results in that case did not depend on Nf 
nor they depend on a magnitude of the chiral condensate in deconfined phase or a value of the quark's mass if it 
would be nonzero. The same remark also applies for analysis of dense matter with p ^ discussed in the previous 
section as long as k = Nf/N — > 0. 

In this section we want to speculate what happens when k = Nf/N ~ 1 by considering hot matter T^O, p ~ 0. 
In contrast with previous analysis we anticipate a very strong dependence from all (previously unessential) parameters 
such as quark's mass m q , number of flavors Nf = kN, magnitude of the chiral condensate etc. We start by 

considering variation of transition properties on quark's masses. For simplicity, consider the limit when all quarks 
have the same and sufficiently large mass m q 3> Aqcd- It is obvious that the first order phase transition anticipated 
in pure gluodynamics will not be effected by presence of Nf sufficiently heavy fermions as they essentially decoupled 
from the system in the limit m q ^> Aqcd- When m q is getting smaller but still sufficiently large m q > Aqcd one 
can easily demonstrate that the structure of 7(2"') remains the same, but the corresponding critical temperature will 
be slowly decreasing with m q as follows 

T c (« 9 6 ) m,) = T c (« = 0)fl-«l. (— ) 4 + -)> rn q ^A QCD , «=^~1 (15) 

\ 11 75 m q 11 735 m q J N 

where the few first coefficients of the expansion 1/ (m q p) k in the instanton background have been explicitly calculated 
long ago [37j . see also recent paper [38| . For our estimates (|15l) we replaced p — > (7rT c ) _1 as a typical value of p where 
the integral J dp converges. The expansion (|15|) can be trusted starting from {m q p) ~ m g /(7rT c ) > 1. The fermion 
contribution is a sub leading effect ~ 1/N; it becomes of order one when k = -jf ~ 1, as expected. 

Now, we want to demonstrate a strong dependence of the transition as a function of the chiral condensate (iptp) 
magnitude. If the chiral condensate does not vanish identically in close vicinity of the phase transition, T > T c as 
holographic model of QCD suggests [l|, one can repeat the corresponding calculations with the following result: the 
structure of 7(T) function as defined in ((TJ remains the same while its coefficients would now depend on dimensionlcss 
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parameters re and the value of the chiral condensate 5 . The numerical values of the critical temperature T c and 
coefficient a would change, however the sharp changes of 9 dependence which is a consequence of a generic structure 
of 7(T) remain the same. Therefore, we expect the first order phase transition to hold in this case in complete analogy 
to previously considered case re — * 0. 

However, we think it is very unlikely for the chiral condensate to remain finite at T > T c when re = Nf/N ~ 1. It 
is much more likely that the chiral condensate vanishes at T > T c when re ~ 1. In this case our analysis based on 
the dilute instanton approximation (jHJ) will be obscured due to the long range interactions between instantons and 
anti-instantons induced by massless quarks, [39|, |40(. This induced interaction becomes crucial even when instantons 
are still far away from each other, and the instanton gas is still dilute. The corresponding estimations for T c and 
studying the properties of the transition in the case re ~ 1 become very model dependent analysis, and we shall 
not elaborate on this issue in the present paper. We anticipate that the transition properties will be very sensitive 
to the quark's masses as the instanton interactions drastically depend on the quark's features in this case. Such 
a sensitivity is consistent with the lattice results which suggest that for vanishing quark masses there will be first 
order phase transition while for physical masses it becomes a smooth crossover, see e.g. recent reviews (UH^]. We 
should emphasize here that our basic principle which relates sharp changes in 9 dependence and transition properties 
still holds for re ~ 1 case. This principle is a simple reflection of the fact that the point where the confinement sets 
in corresponds to the regime where the instanton density suddenly blows up and the instantons dissociate into the 
instanton quarks as mentioned in chapter II and discussed in a more details in (lj and references therein. The large 
number of fermions when re ~ 1 obscures a simple analysis when the critical point can be estimated by approaching 
from deconfined phase where the instanton density is parametrically suppressed ([3]) and the system remains under 
theoretical control up to a close vicinity of T c . In the case of re ~ 1 the corresponding analysis becomes much more 
involved due to the reasons mentioned above. 

Therefore, the main lesson from estimates presented above is as follows. In the case when Nf ~ N we observe a 
great sensitivity of the transition properties on specific details of the system such as quark' s masses, magnitude of 
the chiral condensate, value of re. It is very difficult to make any solid predictions in this situation as they would 
largely depend on underlying assumptions. Such a sensitivity of the transition to quark's properties at re ~ I is in a 
huge contrast with our previous estimates when there is unambiguous prediction for the first order phase transition 
at T ~ I and small chemical potential (|2I3I5P and very large \i ~ \f~N and small temperature (|12I13I14|) when re <C I 
and when all specific quark's details are irrelevant as their contribution is suppressed at least by factor 1 / \f~N . 



V. SUMMARY. FUTURE DIRECTIONS 



We explore the consequences of the assumption that in the large N QCD at Nf -C N confinemcnt-dcconfincment 
phase transition takes place exactly at the point where the dilute instanton calculation breaks down, and therefore 
where 9 dependence must drastically change. This conjecture for T ^ is supported by lattice computations and 
holographic arguments. At large chemical potential we do not have such independent support. However, the basic 
governing principle remains the same, and therefore, our results (|12I13I14|) can be considered as a prediction. The 
most important consequence of our conjecture is observation that the critical \x c is very large, \x c ~ \f~N~ which is 
consistent with fundamentally different arguments presented in ref. [2|. Another important observation is the fact 
that the first order phase transition at re <C 1 holds all the way down from T c ~ \, /i = to fj, c ~ V~N, T = 
as a consequence of the same nature of the transition. In different words, the 9 behavior experience sharp changes 
whenever the phase transition line is crossed. This feature is very robust consequence of our conjecture, not sensitive 
to details of quark's properties such as masses, chiral condensation etc, as they may influence the sub-leading terms 
only. Situation becomes drastically different when Nf ~ N in which case everything becomes very sensitive to details: 
quark masses, chiral condensate, precise value of re, etc. 

A general comment on this proposal can be formulated as follows. Our conjecture which relates two apparently un- 
related phenomena (phase transition vs sharp changes in 9 behavior) implicitly implies that topological configurations 
which are linked to 9 must play a crucial role in the dynamics of the phase transition. For T > T c such configurations 
are well-known: they are dilute instantons with density ~ e~ lN cos 9. We presented arguments in[l[ (see also earlier 
references therein) suggesting that at T < T c the instantons do not disappear from the system, but rather dissociate 



Indeed, in our previous analysis the chiral condensate enters the instanton density as follows ~ (ipip) N f ~ e W'(«ln |{V>^)l) a For re = 

-J- — + this term obviously leads to a sub leading effects 1/JV in comparison with the main terms in the exponent ((TJ- For re ~ 1 this 
terms becomes the same order of magnitude as other contributions in JTJ- 
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into fractionally charged constituents, the so-called instanton quarks. In this sense the phase transition can be under- 
stood as a phase transition between molecular phase (deconfined) and plasma phase (confined) of these fractionally 
charged constituents. The same arguments still hold for the entire phase transition line in (T, /z) plane. A similar 
conclusion on sharp changes in 6 behavior at T = T c was also observed in ref.[24| where the authors studied the D2 
branes in confined and deconfined phases at T ^ 0. The topological objects (sensitive to 9) were identified in ref.[24j 
as magnetic strings. 

Our final remark here is as follows. If the picture advocated in the present work about the nature of the transition 
turns out to be correct, it would strongly suggest that fractionally charged constituents (which carry the magnetic 
charges as discussed in[l| and references therein) may play a very important role in dynamics in deconfined phase in 
close vicinity of the transition < (T — T c ) < 1/N. For large N this region shrinks to a point, however for finite N it 
could be an extended region in temperatures. In this region the instantons are not formed yet, and our semiclassical 
analysis is not justified yet as eq. suggests. However, the constituents in this region are already not in condensed 
form. Therefore they may become an important magnetic de gre es of freedom which may contribute to the equation 
of state, similar to analysis on wrapped monopoles in ref . [43L 13] . The role of these fractional magnetic constituents 
could be even more profound if Nf ~ N where smooth crossover likely to take place jZJ . In this case the region 
of interests is order of one (T — T c ) ~ 1 in Aqcd units in contrast with a narrow region < (T — T c ) < 1/N if the 
first order phase transition takes place. The region above T c is also very interesting from phenomenological viewpoint 
as reviewed in [45J . 
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